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1. Why NNLO calculations are important. 

In the past decade, QCD has become a quantitative science and comparisons of 
experimental data with NLO QCD have become the standard. However, there are 
many reasons why extending perturbative calculations to NNLO is vital in reducing 
the theoretical uncertainty. In the following we list six of them. 


1.1 Renormalisation scale uncertainty 

In many cases, the uncertainty from the pdf’s and from the choice of the renormali¬ 
sation scale give uncertainties that are as big as or bigger than the experimental 
errors. Of course, the theoretical prediction should be independent of hr. However, 
a scale dependence is introduced by truncating the perturbative series. The change 
due to varying the scale is formally higher order. If an observable Obs is known to 
order , 

N 

Obs = '^An{fj,R)a'^{fj,R), 


then. 


-Obs = O (q,''+') . 


c)ln(,,y V * / 

Often the uncertainty due to uncalculated higher orders is estimated by varying the 
renormalisation scale upwards and downwards by a factor of two around a typical 
hard scale in the process. However, the variation only produces copies of the lower 
order terms. 


Obs = AittsifiR) 

+ ^^2 + boAi In ^ ■ 

A2 will generally contain infrared logarithms and constants that are not present in 
Ai and therefore cannot he predicted by varying pR. For example, Ai may contain 
infrared logarithms L up to L^, while A2 would contain these logarithms up to 
fiR variation is only an estimate of higher order terms and a large variation probably 
means that predietable higher order terms are large. 

To illustrate the improvement in scale uncertainty that may occur at NNLO, let 
us consider the production of a central jet in pp collisions. The renormalisation scale 
dependence is entirely predictable, 
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Figure 1: Single jet inclusive distribution at Et = 100 GeV and 0.1 < \ri\ < 0.7 at 
^/s = 1800 GeV at LO (green), NLO (blue) and NNLO (red). The solid and dashed red 
lines how the NNLO prediction if A4 = 0, ^44 = ±^ 1/^2 respectively. The same pdf’s and 
as are used throughout. 

with L = \og{^R/E t). A 2 and ^3 are the known LO and NLO coefficients while 
A4 is the presently nnknown NNLO term. Inspection of Fig. 1 shows that the scale 
dependence is systematically rednced by increasing the nnmber of terms in the per- 
tnrbative expansion. At NLO, there is always a tnrning point where the prediction 
is insensitive to small changes in If this occurs at a scale far from the typically 
chosen values of /i/j, the iT-factor (dehned as iT = 1 + as{iJiB)A^/A2) will be large. 
At NNLO the scale dependence is clearly signihcantly reduced, although a more 
quantitative statement requires knowledge of A4. 

1.2 Factorisation scale dependence 

Similar qualitative arguments can be applied to the factorisation scale inherent in per¬ 
turbative predictions for quantities with initial state hadrons. Including the NNLO 
contribution reduces the uncertainty due to the truncation of the perturbative series. 

1.3 Jet algorithms 

There is also a mismatch between the number of hadrons and the number of partons 
in the event. At LO each parton has to model a jet and there is no sensitivity to the 
size of the jet cone. At NLO two partons can combine to make a jet giving sensitivity 
to the shape and size of the jet cone. Perturbation theory is starting to reconstruct 
the parton shower within the jet. This is further improved at NNLO where up to 
three partons can form a single jet, or alternatively two of the jets may be formed 
by two partons. This should lead to a better matching of the jet algorithm between 
theory and experiment. 

1.4 Transverse momentum of the incoming partons 

At LO, the incoming particles have no transverse momentum with respect to the 
beam so that the hnal state is produced at rest in the transverse plane. At NLO, 
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Figure 2: The average value of (1 — T) given by Eq. 1.1 showing the NLO prediction 
(dashed red), the NLO prediction with power correction of A = 1 GeV (solid blue) and an 
NNLO estimate with ^3 = 3 and a power correction of A = 0.5 GeV (green dots). 

single hard radiation off one of the incoming particles gives the hnal state a transverse 
momentum kick even if no additional jet is observed. In some cases, this is insuf- 
hcient to describe the transverse momentum distribution observed in the data and 
one appeals to the intrinsic transverse motion of the partons conhned in the proton 
to provide an enhancement. However, at NNLO, double radiation from one parti¬ 
cle or single radiation off each incoming particle gives more complicated transverse 
momentum to the hnal state and should provide a better, and more theoretically 
motivated, description of the data. 

1.5 Power corrections 

Current comparisons of NLO predictions with experimental data generally reveal the 
need for power corrections. For example, in electron-positron annihilation, the exper¬ 
imentally measured average value of 1-Thrust lies well above the NLO predictions. 
The difference can be accounted for by a 1/Q power correction. While the form 
of the power correction can be theoretically motivated, the magnitude is generally 
extracted from data and, to some extent, can be attributed to uncalculated higher 
orders. Including the NNLO may therefore reduce the size of the phenomenological 
power correction needed to ht the data. 

Before the calculation of the NNLO contribution it is not possible to make a more 
quantitative statement. However to illustrate the qualitative point, let us take the 
simple example of an observable like (1 — T) which can be modelled by the simplihed 
series, 

(1 — T) = 0.33as((5) + l.Oas(Q)^ -F A^asiQ)^ + (1-1) 

with as{Q) ~ 67 r/ 23 / log(Q/A) and A = 200 MeV. Fig 2 shows the NLO perturbative 
prediction A 3 = 0, A = 0 as well as the NLO prediction combined with a power 
correction, A 3 = 0, A = 1 GeV which can be taken to model the data. If the 
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Kinematic scales Tree One-loop Two-loop 
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Table 1: Ratios of kinematic scales for the scattering of unlike massless quarks, qq q'q’ 

NNLO coefficient turns out to be positive (which is by no means guaranteed), then 
the size of the power correction needed to describe the data will be reduced. For 
example, if we estimate the NNLO coefficient as = 3, which is large but perhaps 
not unreasonable, then the NLO prediction plus power correction can almost exactly 
be reproduced with a power correction of the same form, but A = 0.6 GeV. We are 
effectively trading a contribution of 0{1/Q) for a contribution of 1/log^((5/A). At 
present the data is insufficient to distinguish between these two functional forms. 

1.6 The shape of the prediction 

It is an oft repeated statement that NLO predictions affect the normalisation of the 
leading order prediction but not the shape of the distribution. This is wrong for two 
reasons. First, as indicated in Table. 1, the loop contribution gets new kinematic 
contributions. Of course, each ratio of kinematic scales is further multiplied by 
logarithms and polylogarithms of ratios of scales. Second, as the number of hnal 
state particles increases, the allowed phase space is enlarged - leading, for example, 
to lower values of thrust in electron-positron annihilation and jet production at higher 
rapidities in hadron-hadron collisions. These effects simply cannot be modelled by 
scaling the lowest order cross section. 

1.7 Parton densities at NNLO 

Consistent NNLO predictions for processes involving hadrons in the initial state 
require not only the NNLO hard scattering cross sections, but also parton distribution 
functions which are accurate to this order. 

The evolution of parton distributions is governed at NNLO by the three-loop 
splitting functions, which are not fully known at present, although the expectation 
is that they will be determined shortly, see Ref. [1] and references theirin. 

The determination of NNLO parton distributions requires a global £t to the avail¬ 
able data on a number of hard scattering observables, with all observables computed 
consistently at NNLO. At present, the NNLO coefficient functions are available for 
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the inclusive Drell-Yan process [2, 3] and for deep inelastic structure functions [4], 
These two observables are by themselves insufficient to fully constrain all parton 
species. In particular large transverse energy jet production in hadron-hadron col- 
lisons provides a vital constraint on the gluon at large x and therefore requires the 
determination of the single jet inclusive cross section at NNLO. 

2. Recent progress in the field 

There has been enormous progress in calculating two-loop matrix elements since the 
previous “Loops and Legs” workshop in April 2000. To a large extent this has been 
due to breakthoughs in the calculation of two-loop master integrals with massless 
propagators. In tour de force calculations through the summer of 1999, Smirnov [5] 
and Tausk [6] computed the planar and non-planar double box master integrals with 
four on-shell legs using Mellin Barnes integrals. The other four four-point master 
integrals needed for the on-shell case were much easier to calculate [7, 8, 9, 10, 11]. 

Smirnov has applied the same Mellin Barnes technique to the case where one leg 
is off-shell [12, 13], however, the main break through here is due to Gehrmann and 
Remiddi [14, 15] who have provided expansions for all of the relevant master inte¬ 
grals by making use of differential equations for two-loop four point functions [16]. 
For integrals with singularities through to 1/e^ it proved necessary to introduce 
two-dimensional harmonic polylogarithms in addition to the standard Nielsen poly¬ 
logarithms. Numerical evaluations of the master integrals [17, 18] and analytic con¬ 
tinuations have been provided in Ref. [19]. 

The next major step is to reduce the tensor integrals appearing in Feynman 
diagram calculations to the few master integrals for which series expansions are 
known. The main tool here is integration-by-parts (IBP) [20, 21] which, for two-loop 
four point integrals, generates 10 equations relating integrals with different powers of 
propagators. To this can be added the three Lorentz Invariance (LI) identities [16]. 
Together, these equations represent the fact that loop integrals exhibit Poincare 
invariance. Together, the IBP and LI equations form a linear system that can be 
automatically solved with algebraic methods to relate all tensor integrals to the 
master integrals [16, 22]. 

These technological tools have allowed the computation of a whole raft of two- 
loop matrix elements for scattering processes with up to one off-shell leg as detailed 
in Table 2. 

A vital check on the infrared structure of these two-loop matrix elements is given 
by Catani [37]. In Catani’s formalism, the singular (P(l/e^), (P(l/e^), (P(l/e^) struc¬ 
ture can be straightforwardly be determined while the non-logarithmic coefficient 
H 2 multiplying the (P(l/e) can also be extracted. For some time, the origins of this 
formula were lost in transit, however Sterman and Tejeda-Yeomans [38] have now 
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On-shell Process 


Tree Helicity 

X Two-loop amplitudes 


e'''e —(e+e ) 

qq qq{q'q') 

99 

99 99 

99 77 

77 ^ 77 
qq 9l{ll) 


[23] 

|24, 25] 

[26] 


[27] 

[28, 29] 

— 

[30] 

[33] 

[31, 32] 


Off-shell Process 

e+e- ^ qqg \M] [35, 36] 


Table 2: Current status of available two-loop amplitudes for scattering processes with all 
legs on-shell and with one leg off-shell 

rederived it from the basis of factorisation of inter-jet and intra-jet radiation and 
given a more concrete explanation of how H 2 is produced. 

3. What remains to be done 

Of course knowledge of the two-loop amplitudes is only one of the ingredients needed 
in the construction of a NNLO parton level Monte Carlo. In addition we also need; 

• square of the one-loop 2 —2 amplitudes, 

• interference of tree and one-loop 2 —3 amplitudes, 

• square of the tree-level 2 —4 amplitudes. 

The tree level 2 —4 amplitudes have been known for some time, while the hve point 
amplitudes with zero [39, 40, 41] and one [42, 43, 44, 45] off-shell legs have also been 
worked out. 

These latter two processes contribute when either one or two of the partons are 
unresolved and there is a much more sophisticated infrared cancellation between 

• n and n + 1 particle contributions when one particle is unresolved, 

• n and n + 2 particle contributions when two particles are unresolved. 

At present, a detailed method for subtracting the singularities and numerically eval¬ 
uating the hnite remainders has not been worked out. 

Nevertheless, the prognosis is favourable. If the rate of technical development 
maintains its current pace, it is extremely likely that by the time of the next work¬ 
shop, NNLO Monte Carlo’s will be available for some of the most basic scattering 
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processes, leading to a more complete and quantitative understanding of hard inter¬ 
actions. 
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